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A METHOD FOR CALCULATING THE AERODYNAMIC LOADING 
ON WING— TIP -TANK COMBINATIONS 
IN SUBSONIC FLOW 

By Samuel W. Robinson, Jr., and Martin Zlotnick 


SUMMARY 


An analytical method for calculating the aerodynamic loading on 
■wing— tip-tank combinations in subsonic flow is developed by using a 
simple horseshoe vortex- image system for the case of the tank axis in the 
plane of the wing. 

An illustrative example is given in the appendix, in which wing 
and tip-tank loadings are calculated for three configurations. The 
calculated results are shown to be in good agreement with experimental 
data. 


INTRODUCTION 


With external wing tip tanks coming into more general usage, an 
analytical method for calculating the effects of mutual interference 
between wing and tip tank will be valuable. This interference may 
cause large differences in wing torsion and bending moment from the 
wing-alone values and may introduce significant changes in the perform- 
ance and stability characteristics of the aircraft. 

Reference 1 presents a method for calculating the lift induced on 
the fuselage of a wing -fuselage combination for incompressible flow 
using as a basis for its development the mathematical model consisting 
of a set of horseshoe vortices and their images to represent a wing of 
arbitrary plan form in the presence of an infinite cylinder. This 
method affords a means of determining the distribution of the induced 
lift on the body and of the additional wing lift and is valid for con- 
figurations of varying plan forms and wing-span-— body-diameter ratios. 
In using this method to compute loadings at high subsonic speeds, the 
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Prandtl-Glauert correction for compressibility may be applied in the 
same manner as for a wing alone. 

In this paper, an extension of the work of reference 1 to the 
unsymmetrical case is made to provide a method of calculating the 
desired loadings on the wing— tip-tank combinations in which the tank 
axis lies in the plane of the wing. 

SYMBOLS 

A aspect ratio, b^/s 

a radius of an infinite circular cylinder representing the tip 

tank 

b wing span 

lift coefficient, L/qS 
Cl.j. tank lift coefficient, Lt/qS 

Cj^ tank pitching-moment coefficient, M /gS d^ v 

cj section lift coefficient, Z/qe 

c wing chord 

d diameter of tip tank 

e length of tip tank 

F, f downwash functions for trailing vortex 

G, g downwash functions for bound vortex 

l section lift 

L -cotal lift 

M total pitching moment 

p local pressure 

f 

q free-stream dynamic pressure 
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r^ radius of tip tank 

s semispan of reference horseshoe vortex, b/l6 

S wing area 

U longitudinal component of free-stream. velocity 

u longitudinal component of local velocity 

V volume 

w vertical component of local velocity (downwash) 

x,y,z longitudinal, lateral, .and vertical coordinates, respectively 

x longitudinal ordinate of center of pressure on tank 

AyU) width of an image vortex 

a angle of attack 

F vortex strength 

a correction factor for finite "body length 

Tj correction factor for semi-infinite body 

X taper ratio 

A sweep angle 

Superscripts: 


* 

dimensionless 

with respect to 

b/2 

! 

dimensionless 

with respect to 

s 

T 1 

dimensionless 

with respect to 

r t 

1 f T 

dimensionless 

with respect to 

^max 

(i) 

function of image vortex 
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Subscripts: 

n running index, lift points or horseshoe vortices 

m number of horseshoe vortices representing wing 

ot center line of tip tank 

p index of downwash stations 

g geometric 

t pertaining to tip tank 

w pertaining to wing 

wo pertaining to wing alone 

LE leading edge of wing tip 

T pertaining to wing tip 

u upper surface of body 

l lower surface of body 

wt pertaining to wing— tip-tank combination 

B pertaining to body-alone component 

( ) a slope per radian, d/da 

BASIC CONSIDERATIONS 


For analytical treatment, the wing— t ip-tank combination in subsonic 
flow may be represented by a horseshoe vortex system similar to that 
used in reference 1 for wing-fuselage combinations. The wing is replaced 
by a set of horseshoe vortices and the tip tank, represented by an 
infinite cylinder, is replaced by the images of the wing vortices. 

A mathematically rigorous solution to this idealized problem would 
require calculation of an additional potential within the cylinder, 
since the image vortices alone satisfy the boundary conditions only in 
particular locations on the cylinder. However, as in reference 1, this 
additional potential can be shown to produce zero net lift and moment 
on the cylinder. Therefore, using only the vortex- image system, a 
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straightforward method may he developed for calculating the lift induced 
on a cylindrical body attached to the wing tip. With modifications to 
the vortex- ima ge system, the spanwise lift distribution on the wing 
with tank attached can be calculated approximately. 

In addition to the induced lift component of the loading on the 
tip tank, which is independent of the angle of attack of the tip tank 
itself and which depends only on the spanwise loading of the wing, 
there is a relatively small component of loading due to the angle of 
attack of the body (referred to hereinafter as the "body alone" lift) 
which can be calculated approximately using Munk’s airship theory. To 
treat a practical configuration, it is necessary to correct the induced 
lift on the tank slightly for finite body length, and the body-alone 
lift for the particular tank location. 

The downwash behind the wing of a wing— tip-tank combination can 
also be calculated if desired by using the method of reference 1 as 
modified in this paper . 


CALCULATION OF LIFT AND MOMENT ON TIP TANKS 


In this section, the derivation of the expression for the induced 
lift on an infinite cylinder due to an external horseshoe vortex and 
its image within the cylinder is presented. The lift of each of a 
group of such vortices is superimposed to calculate the lift induced on 
a tip tank. This is followed by a calculation of the body-alone lift 
on a slender body of revolution, with corrections for the particular 
configuration. Added together, these two components give the total 
lift and moment on the tip tank. 


Induced Lift 

Vortex- image system .- A horseshoe vortex of width yg - y^ outside 
of an inf inite cylinder of radius a is shown in figure 1. An image 

2 „ 2 

horseshoe vortex of width - — lies within the cylinder, satisfying 

y i y 2 

the boundary condition of zero flow normal to the cylinder wall only at 
x = ±oo j in the yz -plane at x = 0, and in the plane z = 0. Although 
in order to satisfy the boundary conditions everywhere else it would 
be necessary to superimpose an additional potential on the system, this 
additional potential has been shown in reference 1 to have zero net 
lift and pitching moment so that it need not be considered in this 
section. 
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Calculation of induced lift on a cylinder due fro a vortex and its 
image . - The induced lift on a section dy of the cylinder in figure 1 
is the integral of the pressure difference on the upper and lower sur- 
faces of the section dy from x = -<» to x = <». 

At two corresponding points on the upper and lower surfaces of the 
cylinder the local pressures are 

Pu " \ p( u + ^u ) 2 


Pi = \ P( u + Au l) 2 


The lift on the section dy is 


£ - (>u - *»)«* 


Now the vortex flow above the cylinder is equal and opposite to the flow 
at corresponding points below the cylinder; therefore, 

Au u ■ ~ Au l 


Expanding and simplifying the expressions for p u and pj 


yields 


dL 

ay 



2pU Au dx 


( 2 ) 


The velocity induced on the cylinder surface by a vortex extending from 
y-]_ to y 2 at (x,y,z), V is given by the Biot-Savart law as 

yi-ya 
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Since 


dL 

dy 



2pU Au dx 


(6) 


substitution of equation (5) for Au into equation (6) and integration 
yields 


dL 

dy 


pur 

it 


tan 


.1 y 2 - y 


tan 


-i y i 


+ tan 


a2 

z 



= pUr — tan"^ Z ( y l ~ 

* a + y l y 2 " y ( y l + y 2) 


(7) 


Equation (7) is the equation for the lateral distribution of lift induced 
on the cylinder. This can be written as 

-A- & . 2 tan -i ik, rg) (7a) 

<1 1 dy * a2 + y l y 2 “ y ( y l + y 2) 


since 



Integrating to get the total lift induced on the cylinder yields 



= pUP Ay^ 1 ) 


( 8 ) 
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Since the longitudinal distribution of Au is symmetrical about 
the plane x = 0, the location of the bound leg of the horseshoe vortex, 
the longitudinal center of pressure of the induced lift must be on the 
same longitudinal ordinate as the bound leg. 

Superposition of lift due to several vortices .- In the more general 
case of the wing— tip-tank combination, the wing is replaced by several 
horseshoe vortices of convenient width, bound legs centered on the 
quarter-chord line of the wing, or the actual center-of -pressure loca- 
tion of the section, if known. The tip tank is represented by an infi- 
nite cylinder containing images of the wing vortices. A typical repre- 
sentation is shown in figure 2. 

The lift induced on the cylinder by the complete wing vortex- image 
system can be calculated by superposing the lifts due to each of the 
vortices and their images . The lift induced on a cylinder by a group 
of m horseshoe vortices representing a wing then is 


L 


m 


pur r 


n=l 


y n - s. 


n 



(9) 


where y n is the center of the nth wing vortex, s n is the semispan 
of that vortex, T n is the strength of that vortex, a is the radius 

of the cylinder, and — — — is Ay^ 1 ), the width of the 

y n " s n y n + s n 

image vortex. In terms of the wing parameters and the dimensionless 
spanwise loading coefficient cjc*, the equation for the lift -curve 
slope of the tank is 


(C U ) = A. f- Ac*) /_*£ a!_J 

V a 2b { H y n ' s n + s n, 


(9a) 


where cic* = — and y„ is measured from the center of the tank. 

4 bU n 

The moment arm of the induced lift of each vortex- image system x 
is measured from the center of moments to the bound leg of that system, 
and the total pitching moment induced on the cylinder is 


M 


m 


m 


= 2Z ^n^n = 21 pUrn 


n=l 


n=l 



— — ) 

y n + s n 1 


*n 


( 10 ) 
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and the moment coefficient is 



(10a) 


(The moment coefficient may he made dimensionless with respect to any 
length; however, in this report, for convenience in comparing calculated 
and experimental data, the tip-tank maximum diameter is used.) 

Correction for finite length of tank .- To the preceding calcula- 
tion of the induced lift, an approximate correction for the finite 
length of the tip tank may he made hy using the same correction as that 
employed in reference 1 for fuselages. The portion of the tank down- 
stream from the wing may he considered an infinite cylinder in all 
cases; whereas the portion ahead of the wing may he treated as a hody 
of the slenderness ratio of the tank forehody. The correction factor 
should he half that for an infinite cylinder (unity) added to half that 
for the forebody (a). The correction factor tj is 


n - |(i + <0 


where a is the correction factor given in reference 1 for ellipsoids 
of revolution. The correction factor T) is plotted as a function of 
slenderness ratio in figure 3* It can he seen that if this factor is 
neglected, the error is less than 5 percent for most practical tip-tank 
slenderness ratios. 


Calculation of Body-Alone Lift and Moment 

Body-alone lift .- The induced lift on the tip tank is independent 
of the angle of attack of the tank and is dependent only on the wing 
lift distribution. There is, however, a component of the total lift on 
the tip tank due to its angle of attack which is referred to as the 
body-alone lift. In calculating this component of the total lift, the 
effect of upwash of the wing on the angle of attack will he neglected 
because the hody lift is small compared to the induced lift. Since the 
wing tends to make the flow parallel to the tip-tank axis in the region 
of the function of the two, there is negligible hody -alone lift on that 
section of the tank. Because of the tip-tank wake, the lift on the sec- 
tion downstream from the location of maximum tank diameter is small, 
and for practical purposes, negligible. Therefore, the concern is 
mainly with the hody lift over the section of the tank ahead of the 
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■wing. An additional component of lift due to viscous effects which 
varies with the square of the angle of attack is not considered, since 
the angles of attack are assumed to be small. 

The body-alone lift is calculated simply from Munk's airship theory 
(ref. 2) and is expressed as 

Lb = !(K 2 - q d LE 2 a (ll) 


where Kg - K-^ is given in reference 2 as a function of the slenderness 
ratio, and djjg is the diameter of the body at the wing-tip leading 
edge. 


Body -alone moment .- The pitching moment for a slender body of revo- 
lution at an angle of attack is given by reference 2 as 


M = 2qaV(f% - 


( 12 ) 


vhere V denotes the volume of the body of revolution and M is taken 

about the centroid of the body- The body -alone moment of the tank vill 

be 

M = 2qaV’ (12a) 

where V' denotes the volume of the portion of the tank ahead of the 
wing-tip leading edge except when the maximum diameter of the tank occurs 

ahead of the wing-tip leading edge. In this case V 1 is the volume of 

the portion of the tank ahead of the tank wHYiimiw diameter. The moment 
reference axis is at the wing -tip leading edge except when the tank 
maximum diameter is used, in which case the reference axis is at the 
tank maximum diameter. This procedure is illustrated in the numerical 
example given in the appendix. 

Total tank lift and moment .- The total tank lift is the sum of the 
induced and body-alone lifts, calculated from equations (9) and (ll) . 

The total tank moment is the sum of the induced-lift and body-alone- 
lift moments, calculated from equations (10) and (12). These sums can 
be added to the wing lift and moment to produce the total lift and 
moment on the wing— tip-tank combination. 
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Spanwise Wing Lift Distribution 

Outline of method .- In a manner similar to that of reference 1, 
the spanwise wing lift distribution is calculated by equating the angle 
of downwash induced at the three-quarter-chord line expressed in terms 
of the vortex strength to the geometric angle of attack at the three- 
quarter-chord line at several points along the span. The resulting 
simultaneous equations are solved for the vortex strengths. The 
vortex-image system can be used to calculate the downwash at the three- 
quarter-chord line (although the boundary conditions on the cylinder 
are not satisfied everywhere ), provided certain corrections are made in 
the downwash functions. The approximate downwash factors calculated 
herein should be sufficiently accurate for most practical purposes) 
however, more accurate values may be used when they become available. 

Downwash in region near bound vortices .- It was pointed out that 
the boundary conditions on the cylinder of the horseshoe vortex- image 
system are satisfied exactly only infinitely far behind the wing, in 
the plane z = 0, and in the plane x = 0. The additional potential 
necessary to satisfy the boundary conditions everywhere else which con- 
tributed nothing to the tank lift and moment, does induce a downwash. 

In the region near the bound vortex, the downwash due the additional 
potential must be considered. As suggested in reference 1, the effect 
of the additional potential will be approximated by a correction to the 
downwash of the bound vortices which will be assumed to take into account 
the greater part of the effect of the additional potential on the down- 
wash at the three-quarter-chord line. This correction is used in the 
illustrative example in the appendix; however, for wing— tip-tank con- 
figurations, the ratio of tip chord to tank radius is usually large 


(say !£>>,) 


so that the effect on the lift distribution will be found 


to be small and the correction may be neglected if the best accuracy is 
not required. If the correction factor is neglected, the tables of 
reference 3 can be used to obtain most of the downwash functions so 
that the time required for making the calculations will be reduced 
somewhat . 


The correction to the downwash of the bound vortices is obtained 
by treating the flow induced by the bound vortices normal to the wing 
three-quarter-chord line as a two-dimensional uniform rectilinear flow 
past a circular cylinder. The downwash velocity due to the bound vor- 

l 2 

tices w is increased then by [- " — ■) w where y Q t - yp is the 


( ^ ) w 

\y 0 t - yp/ 

' A t * 1 


lateral distance from the center of the tank to the point on the wing 
at which the downwash is to be calculated. 
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Replacing the wing by m horseshoe vortices of convenient width, 
the expression for the downwash angle at any point p(x,y) can be 

written 



where ^CjC*^ = — — and F p p is the downwash function of the trailing 

2 

legs of vortex F n located at x Q ,y n and its images on station Xp,y p , 
and G np is the downwash function of the bound leg of this vortex and 
its images on station x^y^. 

From the Biot-Savart law, F np is the sum of terms having the form 

(1^) 

and G n p is the sum of terms having the form 


■^np ~ 


y v - y n ± s n 


1 + 




*p - *n 

y p - ± s n) 2 + (*p - x n) 2 


1 

+ s n 

Xp - 

\j(y p - y n + s n) 2 + (*p - Xn) 2 


£ 

CO 

1 

g 

1 


- B nf * (*p - *nf 



(15) 


where s n is the semispan of the nth horseshoe vortex. If the correc- 
tion to Gttd, ( — \ , is not used, the values F np + Grip can 

\y 0 t - y P ) 

be obtained from the tables of reference 5 • 
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The local angle of attack of the section must be corrected at p 
as in reference 1 for the presence of the cylinder in the free stream. 

The vertical component of the free stream w is affected by the cylinder 
in the same manner as the bound-leg flow; the correction factor (which 

may not be neglected) is again 


%ff = a < 



The final equation for downwash on the pth station on the wing of a 
wing— tip-tank combination is 





(IT) 


Now at all the downwash stations located midway between the trailing 
legs of the wing horseshoe vortices, an equation for the downwash can 
be written and equated to the effective angle of attack The 

resulting m linear equations with m unknowns can be solved simul- 
taneously for the values of ( c i c *) n * 


It may be noted that the Prandtl-Glauert compressibility correction 
may be used if desired in calculating the loading on the wing at higher 
Mach numbers. Only the calculation of the spanwise lift distribution 
on the wing is changed in this case. All the other calculations 
described in this paper will be unaffected. 

The values of ( c l c *) n can used immediately in equations (7a) , 

(9a), and (10a) to get the lift distribution, total lift, and moment 
induced on the tip tank. 


APPLICATION OP METHOD AND COMPARISON 
OF RESULTS WITH EXPERIMENT 


Used in the application of the method of this report are the wing— 
tip-tank configurations of reference 4, for which experimental data on 
tank lifts and moments were available. These wing— tip-tank configura- 
tions are similar to those in general use: aspect ratio, 5.25; taper 
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ratio, 0.39; sweep, 12.7°; ratio of radius of tip tank to tip chord of 
wing, 0.25; slenderness ratios of tanks, 6:1 and 8:1. The comparison 
of the method with the experimental data for this case should indicate 
whether the method will be valid for most wing— tip-tank combinations. 

In the numerical example in the appendix, the wing is replaced by 
several (that is, 10) horseshoe vortices of convenient width, and the 
tip tank by a system of images of the wing vortices, as shown in fig- 
ure 2(b). Since the induced lift on the tip tank depends strongly on 
the lift distribution near the wing tip, smaller wing vortices were used 
at the tip to give better accuracy in that region. Using equation (17) 
for the down-wash, the spanwise distribution of lift is determined for 

the wing, and then by equations (9)j (10) , (ll) and (12), the lift and 
moment of the tip tank are calculated. It is to be noted that if the 
spanwise lift distribution on the wing in the presence of tip tanks is 
known, the computation is reduced by a considerable amount. 

A comparison of experimental results with those obtained from the 
numerical example in the appendix is summarized as follows (x 1 * 1 is 
positive forward of the tank center of gravity in the following tabulation) 





Theory 

Experiment 

Difference 

/based on experimental^ 
\ value / 

(^ I * wo ) a 

• • • 

• • . • • . 

4.32 

4.06 

6.4 percent 

(° L '' t ) a 

(°o. 

4.74 

4.45 

6.5 percent 

' ' Or 

( CL w°) a 

■ ■ ■ • * • • 

0.097 

O.O 96 

1.0 percent 

Jr? 1 r? 

A 

Central • . . . 

0.0185 

0.0172 

7.6 percent 

Forward .... 

0.0192 

0.0175 

9-7 percent 

Rearward .... 

0.0185 

0.0174 

6.3 percent 

x' ' ' = IS 4 

v u/ a 

Central 

Forward 

0.955 

0.289 

0.884 

0.623 

1.2 percent tank length 
4.1 percent tank length 

Rearward . 

2.000 

2.140 

1.8 percent tank length 


and 


Figures 4 to 6 show comparisons of experimental values of 



with theoretical values plotted from values of 


N). “ 4 
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determined by the method of this paper, assuming the zero lift 


points to he those indicated by experiment . Figures 7 and 8 present 
the calculated spanwise lift distributions over the wing, wing— tip-tank 
combination, and over the tip tank. The small difference between the 


calculated and experimental values of 


( CLw °)a 


ndicates that 


the method of calculating wing— tip -tank interference is good. The 
larger differences between calculated and experimental values of (PLyo) 


and 


( CL vt)a 


can be attributed, in part at least, to the simplified 


representation of the wing by only ten vortices, 


Since the induced lift on the tip tanks is a function of the span- 
wise circulation on the wing, the agreement between calculated and 
experimental values of tank lift indicates that the calculated spanwise 
wing lift distribution is valid, although the experimental distribution 
of lift over the wing— tip-tank combination was not available. 


Since the bcdy-alone lift and moment are only small parts of the 
total lift and moment on the tip tank in the presence of the wing, the 
approximations used in simplifying the calculations of the body-alone 
lift and moment are considered to be good enough for the purposes of 
the analysis in this paper. A possible exception will be the case of 
the extreme forward tank, where the body-alone lift is sufficiently far 
forward to produce an appreciable contribution to the tank pitching 
moment . More accurate methods of calculating the body-alone lift and 
moment might be employed to improve the accuracy of the calculation of 
the total tank lift and moment in the exbreme forward case. 


CONCLUSIONS 


A method of analysis has been developed for calculating the loading 
on wing— tip-tank combinations in subsonic flow, taking into account 
wing— tip-tank interference. Good agreement with experimental results 
for the lift on the combination, the lift on the tip tank, and the center 
of pressure on the tip tank were obtained for three typical configura- 
tions. The error in the calculated center -of -pres sure location on the 
tip tank for the case of a tank located with the nose far ahead of the 
wing-tip leading edge was found to be greater than the errors in the 
other calculated results. Since the loading on the tip tank depends 
largely on the spanwise lift distribution on the wing, the calculated 
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spanwise loading would also be expected to be in good agreement with 
spanwise lift distributions obtained by experiment, although the latter 
were not available for comparison . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX 

NUMERICAL EXAMPLE 
Geometry of Configuration 


In order to illustrate the method, calculations will be made for 
the wing— tip-tank configurations used in the tests reported in refer- 
ence h. Experimental model dimensions are wing span, 92.08 inches and 
tank radius, 1 . 9 T inches and the geometric characteristics as given in 


reference K are: 

A.spect ratio 5.25 

Taper ratio 0.393 

Sweep (quarter chord), deg 12-7 

Tip-tank radius, ^ 7 — 0.0^275 

b /2 


Tank slenderness ratios, — - — 

^max 


Central case 6.0 

Rearward case 8.0 

Forward case 8.0 


Location of tank center of gravity, rearward of wing-tip leading 
edse ’ 

Central case 0.80 

Rearward case (aft case of ref. 4) 1.8l 

Forward case -0.195 

The plan form of the wing is shown in figure 2(a) and the various tip- 
tank configurations in figure 2 (c). 


Calculation of Spanvise Lift Distribution on Wing 

Coordinates of lift and downwash stations .- The wing will be repre- 
sented by five horseshoe vortices on each semispan, three inboard vor- 
tices of span i — j and two outboard of span — . The bound legs 

4 2 8 2 
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are centered on the quarter-chord line and the downwash. points are 
located on the three-quarter-chord line, midway between the trailing 
legs. The origin is at the intersection of the quarter-chord line with 
the axis of symmetry. All distances in this section will be made dimen- 
sionless with respect to s, the semispan of the inboard wing horseshoe 
vortices. The center of each bound leg will be called a "lift point." 
(See fig. 2(a).) The coordinates of the lift points are as follows: 


n 

±y n ’ 

V 

i 

1.00 

0.2251*- 

2 

3.00 

.6761 

3 

5.00 

1.1268 

4 

6.50 

1.4648 

5 

7.50 

1.6902 


The coordinates of the downwash points are: 


P 

> 

+1 

V 

1 

1.00 

2.2486 

2 

3-00 

2.3670 

3 

5.00 

2.4844 

4 

6.50 

2.5739 

5 

7-50 

2.6334 


The coordinates of the image lift points are: 


n 

±y n ’ 

V 

5' 

8.1275 

1.6902 

4' 

8.2736 

1.4648 

3' 

8.3038 

1.1268 

2' 

8.3288 

.6761 


The case where n = 2’ denotes a single vortex which represents in 
addition to the successively smaller images between 

and the center of the cylinder (see fig. 2(b)). 
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Calculation, of downwash . - The downwash function for the trailing 
legs F n p is the sum of f n p (the downwash function of the vortex 

located at lift point x rL ^y n (on p)), f n p' (the downwash function 

of the vortex at lift point Xn>-;y n ); and f n p(^) and f n p(^ ' (the 

downwash functions for the trailing legs of the two image vortices (one 
in each tank) ) 


F = f + f 1 + f (l) + f (l)' 

r np ■‘•np T x np ■ L np x np 


and similarly for the hound legs 

®np = ®np + Snp' + Snp^ + ®np^^ 

Making equations (l^) and (15) dimensionless with respect to 8 “ ^ 
yields 


•np 


y P * - y n ’ + s n’ 


'p 


.1 + 


V - x n ' 


^ - V ) 2 * (y p - - jr n > + s n ')i 


y p ' - yn’ - s n' 


V - Xn' 


1 + 


- V ) 2 + (y P * - y n ' - s n' ) : 


(Al) 


and 


g np - 


1 

yp' " yn' + e n' 

> 

1 

> 

J(y P ' - y n ' + %') 2 + (V - *»') 2 

y P ' - y n ' - s n’ 



^yp’ - y n ' - =n ) 2 ♦ (V - V ) 2 


(A2) 
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Since the loading is symmetrical (r(x^,y n ) = only the 

downwash on one half of the wing need be computed. In order to calcu- 
late the downwash function on station 5 due to all vortices of 
strength substitute the values of x and y n r for n = 5, and 

Xp T and y-p 1 for p = 5 (s n T in this case is equal to 0.5) in equa- 
tions (Al) and (A2)j then 




2 . 633*1 - 1-6902 

^2.6334 - 1.6902) 2 + (0.5) 2 


1 

-0.5 


1 


+ 


2.6334 - 1.6902 

>J(2.6334 - 1.6902 ) 2 + (0.5) 2 


*5,5 = T.5342 

Similarly, 

g 5 ^ 5 = 0.9932 
f 5 ^ 5 ' = -0.0050 
g 5 ^ 5 ' = 0.0003 

fg^ 1 ) = -1.4084 

gj. = 0.1660 
0^0 


and. the functions f- and 

0^0 

negligible . 


The correction factor for the 



g 5 ' of the opposite image are 

bound-leg downwash function 
is 


i 
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so 


and 


1 + 


5,5 


G 


r 5,5 





= 1.1652 


= 6.1208 


= 1.1595 


and the total downwash function is 


F=r 


5,5 + 


1 + 



7. VT361 


The calculation of F^ 5 , G]^ 5 , F 3 ^ 5 , G^ 5 , F 2 ^ 5 , and G 2)5 is 

similar. F-, ^5 and do not include downwash functions due to 

images of r-j_, because these are included as parts of the total image 

represented by . Care must be taken that the proper value of s n 

be substituted in computing f np and g np . 


Distribution of lift over wing .- The downwash at point 5 can now 
be written 



i|^0.1002l(c l c*) l - 0.23465 (c 2 c*) 2 
2.50802(c iG *) 4 + 7.4736 i(c z c*)^J 


0. 67227 ( c i c *) 3 - 


and a similar treatment of the other four downwash 


of s imult aneous equations for 


(51 


m 


( c i°*) n < 


points gives the set 
These are as follows 
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3.0350611* - 1.31827T 2 * - 0.240l8? 3 * - 
0.04918^* - o.o4458r 5 * 

-1. 48562 r x * + 4.2493or 2 * - i.ii66i? 3 * - 
0.12065^* - o.o8244r 5 * 

-0.21090?!* - i.20623r 2 * + 4.45794r 3 * - 
0. 81303 r^* - 0. 25910 r 5 * 

-0.14247?!* - 0. 30619 ? 2 * - 2.75001? 3 * + 
8.38920?^* - 2.43846? 5 * 

-0.10021?!* - 0.23465? 2 * - 0.67227? 3 * - 

2.50802?!,* + 7- 47361 ? 5 * 




(A3) 


where 



Wow the geometric angle of attack, corrected to its effective value in 
the presence of the cylindrical tank by equation (l6), is, since the wing 
is untwisted and the tank is not at an angle of incidence = ctg^ , 

ai = 1.00210 

02 = 1.00397 
a 3 


= 1.01496 


2k 
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= 1.03452 
a 5 = 1.16516 

for a unit geometric angle of attack. Equating the effective angle of 
attack to the downwash angle in equations (A3) and solving for ( c l c *) n 

yields the spanwise lift distribution over the wing 

(cjc*^ = 2.06175 
(c z c*) 2 = 1.96321 
(c z c*) 3 = 1.66579 

(CjC*)^ = 1.38555 
(cjc*) = 1.18883 

o 

These values are plotted in figure 7 • The lift-curve slope of the wing 
in the presence of the tip tank is 

(is) 0 ■ Wa ' I 51 ' 4 ' 58 

n=l 

* 


Calculation of Induced Lift on Tip Tanks 

Induced lift .- When the wing lift distribution in the presence of 
the tip tank is known and represented by a set of finite-width horseshoe 
vortices, the lift induced on an attached cylindrical body at the tip 
is, by equation (9a) , 


where ^Ay*^ i ^) n is the width of the image of the nth wing vortex. 
Substituting (cjc*)^ and (Ay*(i)) n yields 
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(C^) = 1.3125 J^l.188) (0.03185) + (1.386) (0.00465) 

(1.666) (0.00288) + (1.963) (0. 00337)j 


= 0.07318 


for one cylindrical tip tank. The correction factor for the finite 
length of the tank for a nose-section slenderness ratio of 5.4:1 is 
0.945 from figure 3* so that the total induced lift coefficient on one 
tip tank for the central and rearward configurations is 


pL*.) ■ (0.9^5) (0.07318) = 0.0692 


For the tank in the forward position the correction is 0-955 for the 
assumed slenderness ratio of 7*5:1 and 


( 0l t)a 


= (0.955) (0.07318) = 0.0699 


Lateral distribution of induced lift .- Equation ( 7 ) which gives 
the lateral distribution of induced lift on the tank due to a single 
horseshoe vortex and its image can he written 



2s 


>' 1 _ v 2 i 1 

n ^ y 


1 + (V + Sn’HV - s n”) 


2y n ’ 'y 


t I 


where the double prime 1 1 denotes dimensionless with respect to r^ , 
the radius of the tip tank, s n T 1 is the semispan of the nth horseshoe 
vortex and y^ 1 1 is the distance from the center of the tank to the 

lift point of the nth vortex. The section lift at a station y r 1 
the tank then is 


on 



26 


NACA RM L53B18 



f 21 (»*»Tu * 


tan 


2s n " 


JT7 


2, 


n 


1 + (y n “ + 3 n ”)(y n ” - e n ") - 2y n ‘ 'y' ’ 



The values of 


table : 



0.5 are calculated in the following 




0.83 


For y' ' 


= 0.5 
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The lifts at some other stations, which have been calculated similarly, 
are tabulated below and are plotted in figure 8. 


y" 

dL/dy 
q b/2 

• 9 

1.0385 

.5 

.8296 

0 

.6270 

-.5 

.4195 

-.9 

.1763 

-1.0 

.0000 


Calculation of body-alone lift .- In the cases of the central and 
rearward tip-tank configurations, the slenderness ratio of the tank sec- 
tion ahead of the wing -tip leading edge is 3 * 92 : 1 , Kg - K^_ is 0.778, 

and is From equation (ll), the body-alone lift is 

b/2 46.04 


Lb - We, | d LE 2 (K2 - K X ) 

and 



_ it A 
2 4 




* 8 <5 - 25) (4l^) 2(0 - 778) 


= 0.01124 
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For the forward tip-tank configuration, 
nose section is 5.4:1 so that Kg - K-^ 


the slenderness ratio of the 

is 0.873, and ^ 

b/2 46.04 


<S) a - l ( 5 - 25 >(iiii) 2(0 - 8T3 > * °- 01318 


Total lift on tip tank .- The total coefficient of lift on the tip 
tank is the sum of the induced and body-alone lift coefficients, which 
in the central and rearward cases is 



= 0.0692 + 0.0112 = 0.0804 


and in the forward tank case is 

( c L t ) = O.O 699 + 0.0132 = O.O 83 I 

The values of are plotted against in figures 4(a), 5(a), 

and 6(a) for comparison with experiment . 


Total Lift on Configuration 

Wing -alone lift .- For comparison with the experimental data, the 
vortex representation of the wing shown in figure 2(a) was used in cal- 
culating the loading on the wing alone. The section lift coefficients 
so obtained are 


n 


1 

2.0196 

2 

1.9093 

3 

1.5850 

4 

1.2427 

5 

.8996 
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and the total wing-alone lift coefficient is 

(°0 - 2 II t( e » c *)n = *' 32 

a n=l 


Wing— tip-tank combination lift .- The total lift coefficient of the 
wing— tip -tank combination is, for the central and rearward cases. 


(° Lwt )a 


= It. 5791 + 2(0.0804) 


4.74 


so that the total increase in lift on the configuration due to adding 
tip tanks is 


A(C L ) a = 0.42 


Calculation of Moments on Tip Tank 

Moment of induced lift - (central tank) .- The moment of the induced 
lift on the tip tank is the sum of the moments of all components of the 
induced lift. The center of pressure of each component is located at 
the same longitudinal station as the corresponding wing section center 
of pressure. These center-of -pressure locations were estimated using 
the curves of reference 5. The distance of each center of pressure from 
the wing-tip leading edge (dimensionless with respect to the tank diam- 
eter) is (x 1 1 * is positive rearward from the wing -tip leading edge) 


n 

v T T 1 

*m 

5 

+0.27 

4 

-.09 

3 

-.44 

2 

-1.09 
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The moments of the individual components about the wing-tip leading 
edge are 


n 

(•»a 

(Ar* (l >)n 

V 

c z c*(A5r^ i ^) n i n '" 

5 


(0.0318) 

(+0.27) 

+0.0102 

4 

(1.386) 

(.0047) 

(-.09) 

-.0006 

3 

(1.666) 

(.0029) 

(-.44) 

-.0021 

2 

(1.963) 

(.0034) 

(-1.09) 

-0073 


From equation 10 (a ), 

(°Mt) 

a n 


so that about the wing-tip leading edge 



0.0003 


The coefficient of induced lift on the tank is ^ = O.O692. The 

center of pressure of induced lift (rearward of wing leading edge) on 
the tank therefore is 


V ' 


= °« Q 9?3 . = 0.004 
0.0692 


The moment coefficient of the induced lift about the tank center of 
gravity in the central case is 



= (0.0692) (0.800 - 0.004) = 0.055 


and in the rearward case 


= 0.0692(1.810 - 0.004) = 0.125 
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In the forward case it is 


<&) 


= o.o 699 (-o .19 - o.oo4) = -o.oi4 


a 


Body-alone moment . - For the central and rearward cases the value 
of V' in equation 12(a) is equal to the volume of the portion of the 
tank which extends ahead of the wing leading edge. From the dimensions 
given in reference 4 for these two cases, this volume is half that of 
an ellipsoid of revolution having a slenderness ratio of 3-92:1 and a 
maximum diameter of 3-86 inches. ^Note that in the coefficient Cj^_ 

the reference length is still the actual value of d^y = 3-9^- inches.^ 
From reference 2 the value of Kg - Kj_ is 0.778. Then 


N) 


JZ± 


= 2 

a ® ^max 


(Kj - Kj) 


2n 3-86 2 (3-86)(3-92) 


2 (11.22) (144) (3-9^) 

= 0.0145 


-(0.778) 


This moment is about the wing -tip leading edge, and to calculate it 
about the tank center of gravity, the center of pressure of the body 
lift must be known. The center of pressure (ahead of the leading edge) 
is 


(*»)■ 

(N). 


0.0145 

0.0112 


= 1.29 


About the tank center of gravity, therefore, on the central tank 


N) 


a 


= 0.0112(1.29 + 0.80) = 0.023 
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and on the rearward tank 



= 0.0112(1.29 + 1.81) = 0.035 


For the forward case the V' in equation 12(a) is the volume of 
the portion of the tank ahead of the tank maximum diameter. From the 
dimensions in reference 4, the tank maxi mum diameter is located at 
1.27 dmay ahead of the wing-tip leading edge (or 1.07 d^ Y ahead of 

the tank c . g . ) , and the shape of the tank ahead of the maximum diameter 
is that of an ellipsoid of revolution with a slenderness ratio of 5.4:1. 
The value of K % - K^_ in this case is 0.873* Then 



= 2 


2* ^94^ 2 (3*94)(5.4) 

(11.22) (144) (3 *94) 


( 0 . 873 ) 


= 0.0238 


The center of pressure of the body lift is 1.80 d max ahead of the tank 

x ' " . = 1.80. 

coefficient about the tank center of gravity is then 


maximum diameter since 


The pitching-moment 


K). ■ NU 1 - 80 + 1 - fl7) 

= (0.01318) (1.80 + 1.07) = 0.038 

Total tank moment .- The pitching-moment coefficient about the tank 
center of gravity in the central case is (puSj = 0.055 + 0.023 = O.O78, 

in the rearward case is = 0.125 + 0.035 = 0.160, and in the for- 
ward case is = -0 .014 + 0 .038 = 0 .024 . The values of are 

plotted against Ct in figures 4(b), 5(5), and 6(b) for comparison 

■^wo 

with experiment. 
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® Lift point 



(a) Diagram of complete vortex-image system 
for the wing— tip-tank combination. 

Figure 2.- Illustration of wing— tip-tank configuration and vortex-image 
system used in illustrative example. 
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Figure 2.- Continued. 



Cor red ion factor, 


CD 



J tenderness ratio of forward tip-tank section,-^— 

Umax 


Figure 3-- Correction factor for the induced lift on the tip tank. 
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(Id) Tip-tank pitching -moment 
coefficient . 


Figure 4 .- Comparison of theoretical and experimental tip-tank lift and 
pitching-moment characteristics. Central tank. 



(a) Tip- tank lift coefficient. 


(t>) Tip -tank pitching-moment 
coefficient. 


Figure 5.- Comparison of theoretical and experimental tip-tank lift and 
pitching-moment characteristics. Rearward tank. 
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(b) Tip-tank pitching -moment 
coefficient . 


Figure 6.- Comparison of theoretical and experimental tip-tank lift and 
pitching-moment characteristics. Forward tank. 



’l/v/ng with hp tank 
(C L Joe- 4.74 




Wing atone 
(Cl ..Jr, - 4.32 


Dimensionless s panwise or dm are, y* 

Figure 7-- Calculated spanwise lift distributions for a wing with and 
without tip tanks. A = 5.25: A = 12.7°: \ - 0.255: r. * = 0.0428. 





